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Abstract 

The electronic band structure of cubic HfO^ is calculated using an ah initio all-electron self- 
consistent linear augmented plane- wave method, within the framework of the local-density ap- 
proximation and taking into account full-relativistic contributions. From the band structure, the 
carrier effective masses and the complex dielectric function are obtained. The T-isotropic heavy 
and light electron effective masses are shown to be several times heavier than the electron tunneling 
effective mass measured recently. The imaginary part of the complex dielectric function €2(w) is in 
good agreement with experimental data from ultraviolet spectroscopic ellipsometry measurements 
in bulk yttria-stabilized HfO<i as well as with those performed in films deposited with the tetrakis 
diethylamide hafnium precursor for energies smaller than 9.5 eV. 
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Aggressive scaling of the gate lengths and equivalent gate oxide thickness following the 
International Technology Roadmap for Semiconductors (ITRS) 1 is forcing the replacement 
of silicon dioxide as a gate dielectric by high-dielectric constant (high-re) materials to reduce 
leakage currents^, and meet requirements of reliability 3 . Si- related band offsets, permittivity, 
dielectric breakdown strength, interface stability and quality with silicon, and the carrier 
effective masses are key aspects of the new (high-re) oxides that, together with costs of the 
fabrication processes^, have to be considered in extending Moore's law^ to an equivalent 
oxide thickness below 10 nm. 

Cubic hafnium oxide (HfO?) is an important candidate for S1O2 replacement as gate 
dielectric material. It has a dielectric constant £#/o 2 (0) ~ 25 at 300 .fT, which is about six 
times higher than that of silicon dioxide, but more than an order of magnitude smaller than 
that of cubic SrTiO^, £srTio 3 (fy — 300, a competing high-K oxide. However, HfOi has a 
conduction band offset AE c ^Hfo 2 ~ 1-5 — 2.0 eV with respect to silicon, which is more than 
one order of magnitude higher than that of cubic SrTiO^-, AE C) s r Tio 3 ~ O.leV—. The high 
dielectric constant and tunneling barrier (with respect to silicon) of Hf02, together with a 
possible adaptation of the device production line to hafnium oxide thin film growth, make its 
candidate stronger for S1O2 replacement as gate dielectric material. Recently, improvements 
in the growth techniques of thin H fC>2 films lead to a gate dielectric with 0.5 nm equivalent 
oxide thickness 7 . 

A first-principles study of the structural, vibrational, and lattice dielectric properties of 
bulk hafnium oxide in the cubic, tetragonal, and monoclinic phases has been performed 8 , 
as well as a theoretical evaluation of some properties of thin Hf02 films on Si(OOl) 9 . How- 
ever, few results concerning band-structure calculations have been published within the 
non-relativistic approach, or in the relativistic approach including spin-orbit interaction 
effects^ - —. The values of the Hf02 carrier effective masses, which are fundamental to the 
modeling of tunneling currents through Hf02 gate dielectrics^ 1 ^, for example, have not 
been explicitly presented. State of the art full-relativistic calculations of the cubic Hf02 
carrier effective masses and frequency- dependent dielectric function are presented in this 
work. The calculated T-isotropic heavy and light electron effective masses are several times 
heavier than the electron tunneling effective mass measured recently^ 1 ^. The calculated 
frequency behavior of the imaginary part of £Hfo 2 i n the low energy range is in good agree- 
ment with the recent measurements of Lim et al— by far ultraviolet (UV) spectroscopic 
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ellipsometry as well as with those of Edwards^ and of Schaeffer et aM-. 

The ab initio all-electron self-consistent linear augmented plane-wave (FLAPW) method, 
within the local-density functional formalism (LDA) and the generalized gradient approxi- 
mation (GGA), is employed considering both non- and full-relativistic contributions to the 
band structure^. The 5s-, 5p-, 4f-, 5d-, 6s-Hf; 2s- and 2p-0 electrons were treated as 
part of the valence-band states. Cubic fluorite-type HfO"2 belongs to the Fm3m (Oh) space 
(point) group, with Hf in (0,0,0), O in (±1/4, ±1/4, ±1/4), for which the Brillouin zone 
(BZ) is a 14 face polyhedral 30,21 . Through a total energy minimization process within GGA, 
a lattice constant a = 5.16 A was obtained for cubic Hf02- 

Our calculated value a m i n = 5.16 A is close to a = 5.06 A and a = 5.04 A obtained by 
Fiorentini and Gulleri 9 and Demkov^, respectively; a = 5.037 A (LDA) and a = 5.248 A 
(GGA) calculated by Zhao and Vanderbilt 8 . Note that the measured lattice constant a exp 
is in the 5.08-5.30 A range according to the data of Wang et al. 21 . 

Figure 1 presents the band structure of cubic Hf02 along high-symmetry directions in the 
BZ and the total density of states (TDOS), which were calculated both in the non-relativistic 
(dashed lines) and relativistic (solid lines) approaches. The most relevant contributions of 
the relativistic corrections are due to the Hf(p)- and ((i)-derived states, which contribute 
to the low energy pattern (between -10 and -18 eV) of the band structure. A remarkable 
shift, of ~ 6 eV upwards is seen for the if /(p)-derived states due to mass- velocity relativistic 
contributions. The upward shift of the Hf(p) states is followed by a spin-orbit splitting of 
about 2 eV which occurs between —10 and —12 eV. The top of the valence band is mostly of 
0(p) character with some mixture with the Hf(d)- as well as Hf(f)- derived states. HfOi 
exhibits a direct gap at the X-point of 3.65 eV. Demkov^ and Peacock and Robertson^ ob- 
tained 3.4 eV using a plane wave basis set and nonlocal ultrasoft pseudopotentials (CASTEP 
code) within the LDA. Our value for the band gap energy is in good agreement with the 
value of 3.6 eV obtained by Boer and Groot^ via LAPW-LDA calculations within a rela- 
tivistic approach but not including spin-orbit interaction effects. All theoretical values for 
the band gap energy are smaller than experimental, which is due to the well-known under- 
estimation of the energy values of conduction band states in ab initio calculations within 
the local-density functional theory. Through UV ellipsometry spectroscopy measurements 
in bulk yttria-stabilized HfOi crystals Lim et al— have assigned a value of ~ 5.8 eV to 
the 0(p) — Hf(d) band gap. We note that a value of 3.29 eV is obtained from the 
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non-relativistic calculation, and the band gap is of an indirect nature in this case. 

Table I presents the carrier effective masses in the [100] and [111] (T — > L, T — > X, 
X — > r) directions, calculated within and without (see the values in parenthesis) the full- 
relativistic approach. In the X — > W direction, the effective masses values are too high to 
be evaluated within a parabolic fit. In both the valence and conduction bands, the carrier 
effective masses are demonstrated to be highly anisotropic. Relativistic effects are seen 
to be mostly important in the V — L direction for electrons and holes, and in the V — X 
direction in the case of holes. The T-isotropic heavy and light electron effective masses, 
defined as m* so = (8m r _ L + 6mp_ x )/14, have values m* he = 1.350m o and m* e = 0.714m o , 
respectively, where m is the free electron mass. Both m* he and m* le are several times heavier 
than the appropriately named tunnelling effective mass, which was recently estimated as: 
m* = 0.1 mo, through measurements of the temperature dependence of gate leakage current 
and Fowler- Nor dheim tunnelling characteristics in metal/hafnium oxide/silicon structures^; 
and m* = 0.17 mo, through measurements of the direct tunnelling leakage current in n + 
poly-Si gate NMOSFET with Hf0 2 as gate dielectric^. The severe discrepancy between 
the calculated T-isotropic electron effective masses and the experimental tunnelling effective 
mass in Hf0 2 highlights the limitations in taking the latter as a measure of the bulk 
Hf0 2 T-isotropic electron effective mass. The large discrepancy may be related to the 
thin film characteristics and imperfections of the H f0 2 gate dielectrics, as well as the 
influence of Hf0 2 /Si and H f0 2 /metal interfaces on the gate current density in the devices. 
Measurements of the carrier effective masses in cubic Hf0 2 samples have not been performed 
yet, which precludes direct comparison of our theoretical estimates of the carrier effective 
masses with experiment. 

The imaginary part, e 2 (u), of the complex dielectric function was obtained directly from 
full- and non-relativistic FLAPW electronic structure calculations, while the Kramers- 
Kronig relation was used to obtain the real part, £i{uj). They are depicted in Fig. 2 
over the energy range - 16 eV. Relativistic contributions to the real and imaginary parts 
of the dielectric constant are found to be important, being responsible for a considerable 
shift of the main peaks towards higher energies. Moreover, due to the relativistic effects a 
more detailed structure is seen in particular for e 2 (u). Figure 3 shows good agreement of 
the calculated e 2 (u)) with the experimental data reported by Lim et al.— for energies smaller 
than 9.5 eV, and also with those performed by Edwards^ and by Schaeffer et aP,, all using 
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UV ellipsometry spectroscopy. For this comparison the whole calculated spectrum for 62 
was shifted to higher energies, by matching its energy threshold to the experimental value 
of the HjOi gap energy.— 

In conclusion, we have studied the band structure of cubic bulk HjO^ using first- 
principles calculations. Conduction- and valence-band effective masses were obtained and 
shown to be highly anisotropic. Relativistic effects are shown to play an important role, 
reflected in the effective mass values and in the detailed structure of the dielectric function. 
The T-isotropic heavy and light electron effective masses were determined to be several 
times heavier than the electron tunneling effective mass measured recently^^. The calcu- 
lated imaginary part of the dielectric function was shown to agree well with experimental 
measurements for energies smaller than 9.5 eV— ~— . 
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National Research Council (CNPq) under contract NanoSemiMat/CNPq # 550.015/01-9. 
We thank Dr. A. Donegan for a critical reading of the manuscript. 
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TABLE I: Valence- and conduction-band effective masses at relevant symmetry points in the BZ of 
HfOi- Values are in units of the free electron mass, mo- The numbers in parentheses correspond 
to effective-mass values obtained from a non-relativistic calculation. 
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FIG. 1: Band structure of HfC>2 along high-symmetry axis of the cubic BZ. The zero of energy 
(solid horizontal line) was set at the valence band maximum (X-point) for both, the full-relativistic 
and non-relativistic calculations. The total density of states (TDOS) is presented in the righthand 
column of the figure; the main contributions for the DOS peaks are depicted. Full-(non-) relativistic 
calculations are shown by solid (dashed) lines. High-symmetry points of some bands were labeled 
according to their irreducible representations given by group theory. 
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FIG. 2: Real (ei) and imaginary (£2) parts of the HfC>2 complex dielectric function calculated 
within the full-relativistic (solid line) and non-relativistic (dashed line) approaches. The curves 
were plotted using a Lorentzian broadening of 0.1 eV. 
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Fig.3 



FIG. 3: Imaginary (£2) part of the H fC>2 complex dielectric function calculated within the full- 
relativistic (solid line) and non-relativistic (dashed line) approaches. The calculated spectrum for 
62 was shifted to higher energies, by matching its energy threshold to the experimental value of the 
Hf02 gap energy as obtained in Ref. 3]. The experimental data of Lim et al— are depicted by 
full-dots, and those due to Edwards^ and Schaeffer et are depicted by upsidedown full-triangles 
and open-triangles, respectively. 
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